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ABSTRACT 



Context. Supergiant fast X-ray transients are a subclass of high mass X-ray binaries that host a neutron star accreting mass from the 
wind of its OB supergiant companion. They are characterized by an extremely pronounced and rapid variability in X-rays, which still 
lacks an unambiguous interpretation. A number of deep pointed observations with XMM-Newton have been carried out to study the 
quiescent emission of these sources and gain insight into the mechanism that causes their X-ray variability. 

Aims. We continued this study by using three XMM-Newton observations of the two supergiant fast X-ray transient prototypes 
XTEJ1739-302 and IGR J 08408-4503 in quiescence. 

Methods. An in-depth timing and spectral analysis of these data have been carried out. 

Results. We found that the quiescent emission of these sources is characterized by both complex timing and spectral variability, 
with multiple small flares occurring sporadically after periods of lower X-ray emission. Some evidence is found in the XMM-Newton 
spectra of a soft component below ~2 keV, similar to that observed in the two supergiant fast X-ray transients AX J1845. 0-0433 and 
IGR J16207-5129 and in many other high mass X-ray binaries. 

Conclusions. We suggest some possible interpretations of the timing and spectral properties of the quiescent emission of XTEJ1739- 
302 and IGR J08408-4503 in the context of the different theoretical models proposed to interpret the behavior of the supergiant fast 
X-ray transients. 

Key words. X-rays: binaries - stars: individual (XTE J1739-302, IGR J08408-4503) - stars: neutron - X-rays: stars 



1. Introduction 

Supergiant fast X-ray transients are a subclass of supergiant 
X-ray binaries (SGXBs) that host a neutron star (NS) accret- 
ing from the wind of its OB supergiant companion (SFXT, 
Sgueraetal., 2005). In contrast to the previously known su- 
pergiant X-ray bin aries, SGXB s , (i.e., the so-called "classical 
SGXBs", see e.g., IWhite et al.L Il995h . SFXTs are character- 
I* ■ ized by a pronounced transient-like activity. These sources un- 
. £^ | dergo few hours long (as opposed to weeks-months long) out- 
. bursts with peak X-ray luminosities of >10 36 erg/s, and exhibit 
5_i ' a large dynamic rang e in X-ray luminosity (^,10 4 b etween out- 
burst and quiescence; Walter & Zurit a Herasl I2007I) . The origin 
of this extreme variability is still debated, and different models 
have been developed to interpret it. One of these models involves 
a NS accreting matter from the extremely clumpy wind of its 
super g iant companion (lin't Zandl I2005I IWalter & Zurita Herasl 
l2007t iNegueruela et all I2008I) . According to this interpreta- 
tion, the sporadic capture and accretion of these clumps by the 
compact object can produce the observed fast X-ray flares. To 
reach the required X-ray luminosity swing, very high density 
clumps are required (a factor 10 4 -10 5 denser than the homo- 
geneous stellar wind, IWalter & Zurita Herasl l2007l) . Numerical 



simulations of supergiant star winds indicate that these high 
density clumps migh t be pr o duced by instab i lities in the wind 
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_ Oskinova et HI I2007I and refer- 
ences therein). |Bozzo_etal 1 (120081) proposed that the X-ray vari- 
ability of the SFXT sources might be driven by centrifugal 
and/or magnetic "gating" mechanisms that can halt most of the 
accretion flow during quiescence, and only occasionally per- 
mit d irect accretion onto the NS (see also lGrebenev & Sunvaevl 
2007). The properties of these gating mechanisms depend 
mainly on the NS magnetic field and spin period. At odds with 
the extremely clumpy wind model, in the gating scenario a tran- 
sition from the regime in which the accretion is (mostly) inhib- 
ited to that in which virtually all the captured wind material ac- 
cretes onto the NS requires comparatively small variations in the 
stellar wind velocity and/or density, and can easily give rise to a 
very large dynamic range in X-ray luminosity. Yet another model 
was proposed for IGR Jl 1215-5952, so far the only SFXT to 
display regular periodic outbursts, which envisages that the out- 
bursts take place when the NS in its orbit crosse s a high dens ity 
equatorial region in the supergiant's wind dSidoli et al.ll2007l) . 

Except for their peculiar X-ray variability, SFXT sources 
share many properties with the previously known SGXBs. 
Measured orbital periods in SFXTs range from 3 to 
30 days (IGR J 16479-45 14: 3.32 days, IGR J 17544-26 19: 
4.92 days; IGR J 18483-031 1: 18.5 days; SAX J18 18.6-1703: 
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30.0 davs:I.Tain et all l2009t IClark et all l2009t iBird et all [2009; 



IZurita-Heras & Chatvl l2009: Sau eraet all 120071) . and are thus 
similar to those of other SGXBs. The only exception is 
the SFXT IGR Jl 1215 - 5952^ which has an orbital period of 
165 days dSidoli et all 120071) . The rela tively high eccentrici- 
ties inferred for two SFXTs (-0.3-0.7; IZurita-Heras & ChatvL 
2009; Rah oui et all 120081) suggest that these systems might be 
characterized by somewhat more elongated orbits than classi- 
cal SGXBs. The spin period of the NS hosted in these sources 
has been measured only in four cases, the periods ranging 
from 4.7 to 228 s (IGR J1841.0-0536: 4.7 s; IGR J 1843-031 1: 
21 s; IGRJ16465-4507: 228 s; IGR Jl 1215-5952: 18678 s; 
Bam ba et all l2Q0lb ISguera et all 120071 ILutovinov etail [2005; 



Swank et all |2006), and thus being similar to the spin periods 
measured in the classical SGXBs. However, owing to the lim- 
ited duration of most observations, it is not possible to exclude 
that a number of SFXTs have much lo nger spin periods (see e.g., 
ISmith et allll998l;lBozzo et alll2008l) . 

A number of pointed XMM-Newton observations of several 
SFXTs were carried out to study the quiescent emission of these 
sources and gain insight into the mechanism that drives their 
peculiar X-ray activity. During the observation of the SFXT 
IGR J16479-4514, XMM-Newton captured the source undergo- 
ing a very steep luminosity decay from the end of an outburst 
to a much lower state. The latter was interp reted in term s of an 
eclipse of the source b y the companion star (iBozzo et all 12008. 
2009l lJainet all 12009). This observation revealed that in at least 
one case the X-ray variability of a SFXT was due to the obscu- 
ration by the companion star. In the case of IGR J 18483-031 1, 
XMM-Newton help ed identify pulsatio ns in the quiescent X-ray 
flux of this source dGiunta et all 120091) . and thus provided strong 
support for the idea that SFXTs also accrete matter during thei r 
quiescent states (see e.g.. Si doli et all 120071 IBozzo et al .,2008). 

To study the low level emission of SFXT sources, we present 
in this paper quiescent state XMM-Newton observations of the 
prototypical SFXTs XTEJ1739-302 and IGR J08408-4503. In 
Sect. [2l we summarize previous observations of these sources, 
and in Sect. [3] and [4] we present our data analysis and results. 
In particular, we find that the quiescent spectra of these sources 
contain a soft component below —2 keV. We discuss some pos- 
sible interpretations of this component in Sect. [5] A comparison 
is also carried out between the quiescent and outburst spectral 
properties of XTE J 1739-302 and IGRJ08408-4503. Our con- 
clusions are summarized in Sect. [6] 



2. The sources 

2.1. XTE J1 739-302 

XTE J 1739-302 is a SFXT prototype, a nd was discovered with 
RXTE during a bright outburst in 1997 (ISmith et allll998l) . The 
identification of its supergiant companion led to the determi - 
nation of the source distance at 2.7 kpc dRahoui et al l l2008h . 
Several outbursts fr om th is source were de tected later with 
20061). and INTEGRAL (ILutovinov et all 



RXTE (Smith eta 



12005b ISgueraet all I2005L 120061: iBlav et all 120081) . XTE J 1739- 
302 was observed in out burst with Swift /B AT on three occa- 
sio ns, on 2008 April 8 dSidoli et all l2009h . o n 2008 August 
13 dSidoli et alll2009bl) . and on 2009 March 10 dRomano et all 
I2009ir In only the first two cases, Swift slewed to the source and 
observations with the X-ray Telescope, XRT, were carried out. 
During the 2008 April 8 outburst, XRT observed XTE J 1739- 
302 —387 s after the BAT trigger. These data showed that the 
source was rapidly (—1000 s) decreasing in intensity, and the 



X-ray spectrum (0.3-10 keV) could be reproduced well by us- 
ing an absorbed (A r H=13xl0 22 cm~ 2 ) power law (hereafter, 
PL) model (photon index r=1.5 ) The 0.5-100 keV X-ray lu- 
minosity was — 3.0xl0 36 erg/s. Sido li et all d2009l) also per- 
formed an analysis of the Swift broad band (0.3-60 keV) spec- 
trum of XTE J 17 39-302 during this outburst, and found that 
this spectrum could be reasonably well described by using ei- 
ther a power law with a cutoff at high energy (—13 keV), or 
a Comptonizing plasma model (COMPTT in XSPEC). For the 
outburst of 2008 August 13, XRT data were obtained start- 
ing from —390 s after the BAT trigger, and revealed a more 
complex behavior th an that observed during the previous event 
(Sidoli et all l2009bl) . A time-resolved analysis showed that the 
source X-ray spectrum could be fit equally well by using an 
absorbed PL or a black-body (BB) model with constant pho- 
ton index or temperature (T- 1.2, /cTbb^I-8 keV), and a vary- 
ing absorption column density (in the range 3-9 xlO 22 cm -2 ). 
The combined XRT+BAT broad band (0.3-60 keV) spectrum 
could be well fit by using either a model of Comptonization of 
seed photons in a hot plasma (COMPTT in XSPEC) or a BMC 
model. The BMC comprises a BB component and a component 
accounting for the Comptonization of the BB due to thermal 
and/or dynamical (bulk) Comptonization. The 0.1-100 keV X- 
ray luminosity derived from the simultaneous XRT+BAT spec- 
trum was -3.8x 10 36 erg/s. On 2009 March 10, XTEJ1739-302 
again triggered BAT dRomano et all |2009). On this occasion, 
Swift did not perform any quick slew towards the source and 
XRT data were accumulated only ~1.5 h after the BAT trig- 
ger. At this time, the source was already much fainter (X-ray 
luminosity — 7xl0 34 erg/s, 2-10 keV), and the XRT spectrum 
could be reproduced well by using an absorbed power-law model 
(7Y H =4x 10 22 cm" 2 , T=1.2). 

Little is known about the quiescent emission of XTE J 1739- 
302. An ASCA observation in 1999 did not detect the source 
and placed a 3 a upper limit on i ts X-ray lum i nosity of 
8xl0 32 erg/s (exposure time -13 ks. ISakano et all 120021) . A 
—5 ks Chandra observation in 2001 caught the source in a 
relatively low luminosity state (9.7 xlO 33 erg/s) and the X-ray 
spectrum was fit well by using an absorbed powe r-law model 
(N H ~4.2xl0 22 cm" 2 , r=0.62: ISmith et all 120061). B ased on a 
monitoring program with Swift, iRomano et all d2009bl) carried 
out the first study of the long-term variation in the quiescent 
emission from XTE Jl 739-302. They identified three different 
quiescent states of the source, with a 2-10 keV X-ray luminos- 
ity of 10 35 , 2 x 10 34 , and 5 x 10 32 erg/s, respectively. The X-ray 
spectra of these states could be reasonably well described by 
using an absorbed power-law model, with absorption column 
densities and photon indices in the range (0.3-3.6) x 10 22 cm -2 
and 0.5-1.6, respectively. Alternatively, these spectra could also 
be fit by using a BB model with temperatures and radii in the 
range 1.3-1.9 keV and 0.02-0.28 km, respectively. With these 
results at hand, the authors argued that the quiescent emission of 
XTE J 17 39-302 was most likely produced by a spot covering a 
relatively small region of the NS surface, possibly the NS mag- 
netic polar caps. 



2.2. IGRJ08408-4503 

IGRJ08408-4503 was discovered in the Vela region on 2006 
May 1 5 with INTEGRAL during a short flare lasting less than 
1000s(Go tzet al. I l2007l) . Its optical counterpart was later iden- 
tified as the supergiant star HD 74194 located at 3 kpc , thus con- 
firming that this source be longs to the SFXT class dGotz et all 
120071: lMasetti et all I2006T) . IGRJ08408-4503 was observed in 
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XTE .11739-302 (Epic-PN 0BS1) XTE J1739-303 (Epic-PN 0BS3) IGR J05406-4503 (Epic-PN) 
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Fig. 1. Selection of the high background intervals during the two XMM-Newton observations (Epic-PN camera) of XTEJ1739-302 
(OBS1 on the left and OBS2 in the middle) and the observation of IGRJ08408-4503 (on the right). In each case, we reported the 
source lightcurve not corrected for the selection of the good time intervals and not subtracted for the background in the 0.2-12 keV 
(upper panel) and 10-12 keV energy band (middle panel), and the count rate of the total FOV in the 10-12 keV energy band (lower 
panel). In all cases, the time bin is 100 s. Only the observational intervals in which the total FOV count rate in the 10-12 keV energy 
band was below the threshold indicated with a dashed line were considered for the timing and spectral analysis of the sources. 



outburst two times with INTEGRAL, and the combined JEM-X 
and ISGRI spectra were most accurately fit by using a rather flat 
cut-o f f power-law model (T~0, -E cut ~ll-15 keV, Gotz et al., 
120071; iLevder et all 120071) . The absorption column density was 
~0.1 x 10 22 cm~ 2 , comp atible with the inters t ellar v alue in the 
direction of the source (iDickev & Lockman. 1990). The 0.1- 
100 keV luminosity of the two outbursts was 7.6xl0 35 and 
3xl0 36 erg/s. 

IGR J08408-4503 was al so caught in outburst by Swift 4 
times between 2006 and 2009 dRomano et aUl20"08l:ISidoli etaTl 
2009b: iBarthelmv et al., 2009). During the first outburst, which 
occurred on 2006 October 4, Swift/XKY slewed to the source 
~2000 s after the BAT trigger. The combined XRT+BAT spec- 
trum could be reproduced well by using a cut-off power- 
law model (r=0.31, E c ~\\ keV). The column density was 
~0.3xl0 22 cm -2 . The second outburst, which took place on 
2008 July 5, was characterized by complex behavior and com- 
prised two different flares separated by ~40 5 s. The time- 
resolved spectral analysis found that the soft X-ray emission of 
the source could be reasonably well described by using an ab- 
sorbed power-law model with a constant photon index of ~1 
and a variable absorption column density ranging from ^0.5 
to ^llxlO 22 cm~ 2 . The simultaneous XRT+BAT broad band 
(0.3-80 keV) spectrum of the source could be reproduced well 
using an absorbed cut-off power-law model, with E C >14 keV, 
T^l.4, and ^=6.7 xlO 22 cm~ 2 . The third outburst occurred 
on 2008 September 21. On this occasion, Swift /XRT slewed to 
the source ^150 s after the BAT trigger, and time-resolved spec- 
troscopy did not reveal any variation in the absorption column 
density. However, the total XRT spectrum of the observation (ex- 
posure time 1 160 s) clearly revealed a previously undetected soft 
component below ^2 keV. A reasonably good fit to these data 
could be obtained by using an absorbed (7Vh=0.4x 10 22 cm -2 ) 
power-law (T~2.2) component and a BB at the softer ener- 
gies (with a temperature of kT^l.95 keV, and a radius of 
~1.2 km). This BB component was also required to fit the com- 
bined XRT+BAT broad band spectrum of the outburst (a sin- 
gle co mptonization mod el, BMC, only provided a poor fit to the 
data). ISidoli et alj d2009bl) found that the seed photon tempera- 
ture of the BMC component and the BB temperature could not 
be linked to the same value in the fit. This was interpreted as be- 
ing caused by the presence of two distinct photon populations, 
a colder one with a temperature of about 0.3 keV and a radius 



of ^10 km, and a hotter one with a temperature of 1.4-1.8 keV 
and a radius of ~1 km. However, the statistics of the data did 
not allow the authors to distinguish which of these two popula- 
tions was seen directly as a BB and which one ended up being 
seed photons for the thermal Comptonization. The latest out- 
burst from IGRJ08408-4503 was caught by Swift /BAT on 2009 
August 28. Unfortunately, in this case Swift /XRT did not per- 
form any follow-up observation, and thus no detailed spectral or 
timing information was available in the 0.3-10 keV energy band. 

To date, the quiescent emission of IGR J08408-4503 has re- 
mained largely unexplored; the only detection of this source 
during a period of low X-ray activity was obtained by two 
Swift /XRT follow-up observ ations carried out in 2006 May 
22 and 2007 September 29 dLevder et all 120071) . The expo- 
sure time of each of these observations was about ^4.0 ks, 
and only a total of 40 photons could be collected. Assuming a 
Crab-like spectrum, the source X-ray flux was estimated to be 
^2x 10~ 13 erg/cm 2 /s (0.5 - 1 ke V), correspond i ng to a luminos- 
ity of about 2x 10 32 erg/s (Kennea & Campana, 2006). 



3. XMM-Newton data analysis 

For the present study, we used two XMM-Newton observa- 
tions of XTE J '1739-302, and one XMM-Newton observation of 
IGR J08408-4503 . During all three observations, the two sources 
were in quiescence. 

XMM-Newton observation data files (ODFs) were processed 
to produce calibrated event lists using the standard XMM- 
Newton Science Analysis System (v. 9.0). We used the EPPROC 
and EMPROC tasks for the Epic-PN and the two MOS cameras, 
respectively. The event files of the two observations were filtered 
to exclude high background time intervals. The effective expo- 
sure time for each observation and camera is given in Sects. |4~T1 
and 14.21 Source lightcurves and spectra were extracted in the 
0.2-15 keV band for the Epic-PN and 0.2-10 keV for the two 
Epic-MOS camera^. We extracted the background lightcurves 
and spectra from the nearest source-free region to XTE J 1739- 
302 and IGR J08408-4503 . Background and source circles were 
all chosen to lie within the same CCD. The difference in extrac- 
tion areas between source and background was accounted for by 
using the SAS BACKSCALE task for the spectra and the LCMATH 



1 See http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf 
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Fig. 2. XMM-Newton Epic-PN background-subtracted 
lightcurve of XTE J 17 39-302 during the observation car- 
ried out on 2008 October 1 . The upper panel shows the source 
lightcurve in the 0.2-2 keV energy band, whereas the middle 
panel gives the lightcurve in the 2-15 keV energy band (the 
binning time is 300 s). The ratio of the source count rate in the 
two bands, (2-15 keV)/(0.2-2 keV), versus time is shown in the 
lower panel. The time intervals in which no data are plotted 
have been discarded because of high background events. 



task from HEASOFT for the lightcurves. All of the EPIC spec- 
tra were rebinned before fitting so as to have at least 25 counts 
per bin and, at the same time, prevent oversampling the energy 
resolution by more than a factor of three. Given the low count 
rate of the source, the Epic-MOSl and Epic-MOS2 cameras did 
not contribute significantly to the spectral analysis; we therefore 
discuss in this paper only the spectra from the Epic-PN camera. 
Throughout this paper, the errors are given at 90% c.l. (unless 
stated otherwise). 



4. Results 

4.1. XTE J 1739-302 

XMM-Newton observed XTE J 1739-302 twice, on 2008 October 
1 (hereafter OBS1) and on 2009 March 11 (hereafter OBS2). 
The Epic-PN camera was operated in large window mode in the 
first case and in small window mode during the second obser- 
vation. To identify the high background time intervals, we fol- 
lowed the SAS online analysis threacQ and extracted the Epic- 
PN lightcurves for the full field of view (FOV) in the 10-12 keV 
energy band. We discarded time intervals in OBS1 and OBS2 
when the 10-12 keV FOV count rate was higher than 0.3 and 
0.1 cts/s, respectively. This is illustrated in Fig.Q] where we re- 
port the source lightcurves (uncorrected for the background and 
the good time-interval selection) in the energy band 0.2-12 keV 
(upper panel), the contribution of the source during the two ob- 
servations in the energy band 10-12 keV (middle panel), and the 
count rate of the total FOV in the energy band 10-12 keV (lower 
panel). The threshold we imposed in each case is plotted with 
a dashed line in the bottom panel. After the good time inter- 



2 See also http://xmm.esac.esa.int/sas/current/documentation/threads/ 
PN_spectrum_thread.shtml. 
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Fig. 3. The same as Fig. [2] but for the XMM-Newton observation 
carried out on 2009 March 11. The binning time is 300 s. 
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Fig. 4. Hardness-intensity diagram for the two XMM-Newton ob- 
servations of XTE Jl 739-302 (the upper panel is for the observa- 
tion carried out in 2008, the lower panel for the observation in 
2009). The hardness is defined as the ratio of the count rate in 
the hard (2-15 keV) to soft (0.2-2 keV) energy band. The points 
were obtained by using the lightcurves given in Figs. [2] and [3] 
but, wherever necessary, consecutive bins were rebinned so as to 
achieve a S/N>5.5 in the hardness ratio. 



vals were selected, we obtained a total effective exposure time 
of 32 ks and 24 ks for OBS1 and OBS2, respectively. 

The lightcurves of the two observations corrected for the 
good time interval selection and background-subtracted are 
shown in Figs. [2] and [3] In both cases, the source displayed a 
pronounced variability on a timescale of hundreds of seconds, 
with small flares occurring sporadically after periods of lower 
X-ray emission. During these flares, the X-ray flux typically in- 
creased by a factor of ~10-30 (from a few 10 -13 erg/cm 2 /s up to 
~10 -11 erg/cm 2 /s). From the lower panels in Figs.|2]and|3] it is 
apparent that the source emission hardened at higher count rates. 
To investigate this behavior in more detail, we plot in Fig.|4]the 
hardness ratio of the source (i.e. the ratio of the source count rate 
in the hard, 2-15 keV, to the soft, 0.2-2 keV, energy bands) as a 
function of the total intensity (i.e. the total source count rate in 
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Fig. 5. The XMM-Newton spectra of the XTE J 17 39-302 obser- 
vation carried out in 2008. The open circles, open squares, and 
filled circles are from the energy resolved data obtained during 
the time intervals in which the Epic-PN 0.2-15 keV source count 
rate was >0.4, 0.1-0.4, <0.1, respectively. The model used for 
the best fits is an absorbed CUTOFFPL (cutoff energy fixed at 
13 keV). The residuals from this fit are shown in the bottom 
panel. [See the electronic edition of the paper for a color version 
of this figure.] 

the 0.2-15 keV energy band). For these plots, we rebinned the 
lightcurves of Figs. [2] and E] to obtain a S/N>5.5. The scatter in 
the points in Fig. |4]suggests that the relation between the source 
count rate and the hardness is most likely changing from flare 
to flare. However, an overall increase in the hardness with the 
source intensity is apparent from the plots. To test this correla- 
tion, we performe d a linear reg ression that accounts for the errors 
in both variables dPressL 120021) . The slope of the linear function 
determined from this analysis is 1.9±0.4and4.1±0.8forOBSl 
and OBS2, respectively (errors are at 1 a c.l.), thus confirming 
that the correlation is significant (^5 a). To investigate which 
spectral parameters varied with source intensity, we divided the 
lightcurves of the two observations into three intensity levels, 
and accumulated the spectra during the corresponding time in- 
tervals. We selected the intervals in which the 0.2-15 keV Epic- 
PN count rate of XTE Jl 739-302 was < 0.1, 0.1-0.4, and >0.4, 
for both OBS 1 and OBS2. The three spectra extracted from each 

Table 1. Results of the count-rate-resolved spectral analysis 
for the two XMM-Newton observations of XTE J 17 39-302. The 
model used to fit the data is an absorbed CUTOFFPL (we fixed 
the cutoff energy at 13 keV, see text for details). 
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a: in units of 10 22 cm 2 . b: Observed flux in the 0.5-10 keV 
energy band in units of 10~ 13 erg/cm 2 /s. c: Exposure time in ks. 
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Fig. 6. The same as Fig. [5] but for the observation carried 
out in 2009. The model used for the best fits is an absorbed 
CUTOFFPL (cutoff energy fixed at 13 keV). The residuals from 
these fits are shown in the bottom panel. [See the electronic edi- 
tion of the paper for a color version of this figure.] 

XTEJ1739-302 (Epic-PN OBS1) 



0.1 F 




12 5 10 

Energy (keV) 



Fig. 7. The XMM-Newton spectrum of XTE J 17 '39-302 extracted 
by using the total exposure time of OBS 1 . The model used for 
the best fit is an absorbed CUTOFFPL (we did not fix here the 
cutoff energy). The residuals from this fit are shown in the bot- 
tom panel, whereas the middle panel shows the residuals from 
the fit obtained by using an absorbed CUTOFFPL with the cut- 
off energy fixed at 1 3 keV. 



observation were then fitted together with an absorbed cutoff 
power-law model (hereafter CUTOFFP L, see Fig s. [5]and[6]i. We 
fixed the cutoff energy at 13 keV CSid oli et all 120091) . and found 
that this model provided a closer fit to the data than a simple 
absorbed PL model. We report the results of the fit with the 
absorbed CUTOFFPL model in Table Q] (we also attempted to 
fit the data with an absorbed BB model, but this gave a poorer 
result for some of the rate-resolved spectra, with x 2 >1.2). For 
OBS1, the rate-resolved spectral analysis demonstrated that the 
power-law photon index changed from 1.1 ±0.1, in the highest 
source count-rate spectrum, to 1.8±0.3, in the spectrum corre- 
sponding to the lowest source count rate. No significant vari- 
ation in the absorption column density was measured, and the 
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Fig. 8. The same as Fig. [7] but for the observation carried out 
on 2009 March 1 1 . Here the best-fit model to the data was an 
absorbed CUTOFFPL model (the cutoff energy was allowed to 
vary during fitting). 
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Fig. 9. An example of an unfolded spectrum of the XMM-Newton 
observation of XTEJ1739-302 carried out in 2008 (OBS1). Here 
the best fit is obtained by using an absorbed CUTOFFPL plus a 
MKL component (see Sect 14. U . 



best fit value of jVh was consistent (to within the errors) with 
being constant at a level of ~2.6xl0 22 cm~ 2 . A similar anal- 
ysis carried out for OBS2 yielded very similar results. In this 
case, the power-law photon index decreased at the higher source 
count rate from 1.8±0.3 to 1.0±0.1, and the absorption column 
density remained constant (to within the errors) at a value of 
~4.0x 10 22 cm~ 2 . As the measured changes in the spectral pa- 
rameters were not dramatic, we also extracted the X-ray spec- 
tra of the source by using all the data from each observation. 
These spectra are shown in Figs. [7] and [8] A fit to these spec- 
tra with a single absorbed PL or BB model gave a relatively 
poor xled (^ L3 ' d.o.f.=175, 184 for OBS1 and OBS2, respec- 
tively) fit, and a slightly better result was obtained using the same 
CUTOFFPL model discussed before for the rate-resolved spec- 
tra (cutoff energy fixed at 13 keV, ^ red <\2, d.o.f.=175, 184 for 
OBS1 and OBS2, respectively). However, even in this case, the 
reduced \ 2 remained significantly higher than 1 and, especially 



in OBS1, some structures are apparent in the fit residuals (see 
Fig.0. 

By leaving the cutoff energy free to vary in the fit, a sig- 
nificantly closer fit was obtained, with a x 2 ed /d.o.f.=l. 02/174, 
1.09/183 for OBS1 and OBS2, respectively. By using the F- 
test, we found that this improvement was highly significant for 
OBS1 (5.3 a) and somewhat less significant for OBS2 (3.5 a). 
In both cases, the derived cutoff energy turned out to be much 
lower (~4 keV, see Table |2|i than that measured previously 
when the source was in outburst (~13 keV). Such a low value 
for the cutoff energy might not be unlikely for an X-ray pul- 
sar (see Sect. [3). In the case of OBS1, when an additional 
spectral component to that of the CUTOFFPL appears to be 
clearly significant, we also tried to investigate the applicabil- 
ity of other spectral models. We again fixed the cutoff energy 
of the CUTOFFPL component at 13 keV and tried first a phe- 
nomenological model including an additional black-body com- 
ponent (BB) at lower energies (<2 keV). We found that both 
a BB with kT~l-2 keV and a small emitting area (few hun- 
dreds m 2 ) and a much colder (kT^O. 1-0.3 keV) BB with an 
emitting radius of ^100 km could reproduce the spectrum rea- 
sonably well. The interpretation of both of these models for 
the soft X-ray excess faces difficulties: a BB emitting area of 
about one hundred km radius would be much larger than the 
NS, while a small emitting spot on the star surface might also 
be unlikely in the case of a low-luminosity wind-accreting NS 
(see Sect. [5j. We therefore tried to fit the spectrum of OBS1 
using models that provide to alternative physical interpreta- 
tions. We adopted first a model comprising a CUTOFFPL and 
a MKL component, and then tried a second model in which 
we included the effect of partial covering (ZPCFABS in XSPEC) 
on the CUTOFFPL component. We note that these two mod- 
els were suggested to fit the quiescent XMM-Newton spec tra of 
the SFXTs AX J 18 45-0433 (IZurita-Heras & Wa lter. 2009|) and 
IGR J 1 6207-5 1 2S0 dTomsick et al.Ll2009t) . In these cases, the au- 
thors suggested that the MKL component might represent the 
contribution due to the shocks in the stellar wind around the 
NS, whereas the partial covering may be caused by the obscu- 
ration of the NS by clumps in this wind. Finally, we also at- 
tempted to fit the spectrum of OB SI by using the COMPTT 
and the BMC mod els, as suggested by ISidoli et all (120091) and 
ISidoIT et al.ld2009bl) . These two models provided a good fit to the 
data, but in the case of the BMC we found that most of the model 
parameters were only poorly determined in the fit, and thus we 
do not discuss this model further in the case of XTE J '1739-302. 
The results of the fits obtained with all the other models dis- 
cussed in this section are reported in Table [2] We checked that 
a fit to the spectrum of OBS2 with all these models would give 
results comparable to those obtained for OBS 1 . All the spectral 
models reported in Table [2] provided reasonable values of the 
fit parameters and very similar xled ( note mat a H the models 
have the same number of free parameters). In the best case, we 
obtained %^ e(i /d.o.f.=0.98/173, and thus the significance of the 
improvement of these fits, with respect to that obtained with a 
simple absorbed CUTOFFPL model and the energy cutoff fixed 
at 13 keV, was of 5.6 a. As an example, we show in Fig.|9]the un- 
folded spectrum of OBS 1 obtained by using a MKL component 
to fit the low-energy excess. The interpretation of these results is 
discussed in Sect. [5] 

We also performed a Fourier analysis of the lightcurves of 
OBS1 and OBS2, in search of coherent pulsations, by using 



3 This source was clas sified as an intermediate SFXT by 
IWalter & Zurita Herasl d2007h . 



E. Bozzo et al.: XTE J1739-302 and IGR J08408-4503 in quiescence 



7 



Table 2. Spectral fits of the XMM-Newton data of XTEJ1739-302. Below, iV H indicates the absorption column density, T the 
power law photon index, E cut the cutoff energy of the CUTOFFPL component, and Rbb the BB radius. We indicate with kT 
the temperature of the BB component, the temperature of the COMPTT seed photons, or the temperature of the optically thin gas, 
depending on the model. kT c is the temperature of the Comptonizing electron region in the COMPTT model, and r its optical depth. 
We also report the inferred values of Nmkl, which is the normalization of the MKL, Nn 2 , the additional absorption column density 
predicted by the partial covering model, and 0</<l the covering fraction. Finally, we give the source flux in two different energy 
bands (0.5-1.5 keV, 1.5-10 keV), and report also the X-ray flux in the 0.5-10 keV energy band for completeness. 









OBS1 






OBS2 


Model 


CUTOFFPL 


BB+CUTOFFPL 


COMPTT MKL+CUTOFFPL 


ZPCFABS*CUTOFFPL 


CUTOFFPL 






1 7+ ' 4 

— 0.2 

(3.7±0.4) 


1.3±0.2 


3.4±0.4 


1.9±0.4 


7 5 + O.5 
— 0.2 

(3.4±0.2) 
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2+0.3 

(1.3±0.1) 


(1.5±0.1) 


— 


1.5 ±0.1 


1.6±0.1 


l +0 f3 

"•'-0.3 

(1.19±0.07) 


Scut (keV) 


3 5+0.9 
(13 fixed) 


13 (fixed) 


— 


13 (fixed) 


13 (fixed) 


3-8io:e 
(13 fixed) 


kT (keV) 




1.2±0.1 
(0.13±0.01) 


0.8±0.1 


0.15±0.2 
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6 5+ 27 

— 1.6 








A^MKL 








03+1.7 
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3 6+ 15 




/ 










0.7±0.1 




rpd 

r 0.5 -1.5 keV 


1 7+ - 1 

'■'-0.5 

(1.4±0.2) 


i 7 + 0.5 
'-0.7 


i 7 + O.6 
'■'-1.4 


1 7+ 01 

'•'-0.9 


1 s+°- 3 

LS -0.5 


1 1+ - 1 
(0.8±0.1) 


^1.5-10 keV 


25 +o.2 
"-1.1 
(2.5±0.2) 


7 5 + 0.3 
"-1.4 


25+0.5 

"-2.3 


7 5 + 0.3 
"-0.5 


7 s+0.3 
"—0.8 


3 9+04 
(3.9±0.5) 


-^0.5-10 keV 


7 =r + 0.3 
" — 1.0 

(2.5±0.3) 


7 =r + 0.3 
" — 1.4 


25+0.5 

" — 2.0 


2g + 0.2 
"-0.5 


7 5 + 0.3 
" — 0.9 




X r 2 cd /d-o.f. 


1.02/174 
(1.20/175) 


1.02/173 
(0.99/173) 


0.98/173 


1.02/173 


0.98/173 


1.09/183 
(1.16/184) 



a: in units of 10 22 cm 2 . 6: For a distance of 2.7 kpc. c: Absoiption column density of the ZPCFABS component, cl: Observed flux in units of 10 14 erg/cm 2 /s. 
e: Observed flux in units of 10 -12 erg/cm 2 /s. 
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Fig. 10. XMM-Newton Epic-PN background-subtracted 
lightcurve of IGRJ08408-4503 during the observation of 
2007 May 29. The upper and middle panel shows the source 
lightcurve in the 0.2-2 keV and 2-15 keV energy bands, respec- 
tively. The binning time is 300 s. The ratio of the source count 
rate in the two bands versus time is shown in the lower panel. 



the method described in [Israel & Stella (1996). No significant 
(above 3a level) signal was detected in either observation. The 
corresponding 3a cl. upper limits to the pulsed fraction (defined 
as the semi-amplitude of the sinusoid divided by the source av- 
erage count rate), were then comp uted according to the method 
described in Vaugh an et alj (fl99l . In OBS1, upper limits at a 
level of 30%, 20%, and 40% were inferred in the 0. 1-0.2 s, 0.2- 
50 s, and 50-150 s period range, respectively. In OBS2, we de- 
rived upper limits at a level of 20%, 30%, and 35% for periods 
in the range 0.03-20 s, 20-50 s, and 0.02-0.03 s, respectively. 

4.2. IGRJ08408-4503 

XMM-Newton observed IGR J08408-4503 on 2007 May 29, with 
the Epic-PN camera operating in full frame. To identify the high 
background time intervals we followed the same technique de- 
scribed for XTEJ1739-302 (see Sect, and Fig. [Q. We ex- 
tracted the Epic-PN lightcurve for the full field of view (FOV) 
in the 10-12 keV energy band, and set a threshold on the full- 
FOV count rate in this energy band of 0.45 cts/s. The total ef- 
fective exposure time after the good time interval selection for 
IGR J08408-4503 was 26 ks. 

From the lightcurve of the observation (Fig.fTuT). it is appar- 
ent that the variability in the quiescent state of this source was 
rather similar to that of XTEJ1739-302. In particular, the lower 
panel of Fig. [10] shows that the hardness ratio of IGRJ08408- 
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Table 3. The same as Table[2] but for the time-resolved spectra of IGRJ08408-4503. In all cases, the cutoff energy is fixed at 1 1 keV 
(see Sect 14.2b . For the BMC model, we reported below the value T=a+1, where a is the Comptonization efficiency of the model, 
the value of the "illumination parameter" log (A), and the normalization A'bmc- 
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a: For a distance of 3 kpc. b: Observed flux in unit of 10 13 erg/cm 2 /s. c: Exposure ti 
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Fig. 11. Hardness-intensity diagram for IGRJ08408-4503. The 
diagram was created by using the 0.2-2 keV and 2-15 keV 
lightcurves of Fig.[l0] but rebinned so as to achieve a S/N>5.5 
in all the hardness ratio values. The hardness is defined as the 
ratio of the count rate in the hard (2-15 keV) to soft (0.2-2 keV) 
energy band. 



4503 increased with the source count rate. Figure QT| shows 
the hardness-intensity diagram of IGRJ08408-4503, obtained 
with the same technique described in Sect. 14.11 In this case, the 
scatter of the points is somehow less evident than in the case 
of XTE J 1739-302 and a linear fit to the data required a slope 
of 1.9±0.2. To investigate the origin of the variability in the 
hardness ratio of IGRJ08408-4503, we extracted three differ- 
ent spectra during the time intervals of the observation in which 
the source count rate was >0.2, 0.1-0.2, and <0.1 (hereafter 
spectra A, B, C). A fit to these spectra with a simple absorbed 
BB or PL model provided unacceptable results (x;r cd >l. 5-5.0, 
d.o.f. =28-44). A CUTOFFPL model with a fixed E cut =U keV 
(see Sect.[T) provided tighter fits to the three spectra (Xrcd^l-^ - 




1 2 

Energy (keV) 

Fig. 12. The same as Figs. [5] and [6] but for the observation of 
IGR J08408-4503 . Here the open circles, open squares, and filled 
circles refer to the spectrum accumulated during the time inter- 
vals of the observation in which the Epic-PN 0.2-15 keV source 
rate was >0.2, 0.1-0.2, <0.1 counts/s respectively. The model 
used for the best fits comprises a MKL plus a CUTOFFPL com- 
ponent (E cut fixed at 1 1 keV). The residuals from these fits are 
shown in the lower panel. The middle panel shows the residu- 
als of the fits obtained by using a simple absorbed CUTOFFPL 
model (cutoff energy fixed at 1 1 keV). [See the electronic edition 
of the paper for a color version of this figure.] 



than 1, the value of the absorption column density measured 
from the spectra B and C was unreasonably low (compatible 
with zero), and some structures were apparent in the residuals 
from the fits at energies <2 keV (see Fig. [12). A CUTOFFPL 
model with a free to vary E cut improved only the fit to spec- 
trum A (x r ! cd /d.o.f=l. 03/39), whereas the results of the fits to 



1.6, d.o.f=28-44). However, the x^d was still significantly larger the spectra B and C were almost unchanged. To fit the A, B, 
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Fig. 13. The XMM-Newton spectrum of IGRJ08408-4503 ex- 
tracted by using the total exposure time of the observation car- 
ried out on 2007 May 29. The model used for the fit is an ab- 
sorbed CUTOFFPL plus a MKL (the cutoff energy is fixed at 
11 keV). The residuals from the fit are shown in the bottom 
panel. The middle panel shows the residuals of the fit obtained 
by using only an absorbed CUTOFFPL model (cutoff energy is 
fixed at 1 1 keV). 



and C spectra with the same model, and obtain at the same time 
reasonable values of the fit parameters (e.g. an TVh at least com- 
parable with the Galactic value in the direction of the source), 
we assumed tha t E cut =l 1 keV (as it w as found when the source 
was in outburst, iRomano et all l2008h and added a second spec- 
tral component to account for the residuals below ^2 keV. We 
found that a significant improvement in the fits could be obtained 
by adding a MKL or a BB component. A single BMC model also 
provided a reasonable fit to the data. The results of all these fits 
are reported in Table[3] The significance of the improvement ob- 
tained by using these two-component models instead of a simple 
CUTOFFPL is 4, 3, and 4.1 a, respectively, for the spectra A, B, 
andC. 

By looking at the results of the fits in Table [3] we conclude 
that the increase in the source hardness ratio with the count rate 
can be most likely ascribed to a variation in the photon index, 
r, rather than a change in the absorption column density or in 
the properties of the soft component. Indeed, the value of N-r, 
as well as the temperature and normalization of the BB and the 
MKL components, are compatible with being constant (to within 
the errors) in the three spectra. That most of the changes in the 
X-ray flux of the source occurred in the 1.5-10keV energy band, 
where the contribution of the CUTOFFPL component is much 
greater than that of the soft component (see Fig. [9j, added sup- 
port fot the above conclusion. 

To further constrain the soft component more tightly, we also 
extracted the source spectrum by using the total available expo- 
sure time of the XMM-Newton observation. A fit to this spec- 
trum with a simple BB or a PL model did not provide an ac- 
ceptable result (x r ! ed /d.o.f=4.9/101, 1.8/101, respectively for the 
BB and the PL model). A CUTOFFPL model only marginally 
improved the fit (xjr cd /d.o.f=l. 5/100) and some structures re- 
mained present in the residuals (see Fig. fl3l . To obtain an ac- 
ceptable fit to the data, we thus used the same spectral mod- 
els that we adopted for the rate-resolved analysis. All these 
models provided an equivalently good fit to the data. The re- 



Table 4. The same as Table [3] but for the spectrum of 
IGR J08408-4503 obtained by using the total available exposure 
time of the XMM-Newton observation. 
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26 


26 
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a: For a distance of 3 kpc. b: Observed flux in units of 10 13 erg/cm 2 /s. 
c: Exposure time. 



suits of these fits are given in Table [4] and discussed in de- 
tail in the following section. We note that when fitting the to- 
tal spectrum of IGR J08408-4503 we discarded data in the en- 
ergy range 0.4-0.6 keV, as we noted that in this energy range the 
background was rather high. Including these points does not af- 
fect the best-fit values of the model parameters, but indicates 
that the fit with a CUTOFFPL+BB model is slightly prefer- 
able (x r ! ed /dof = 1 - 0/104 ) than 1116 CUTOFFPL+MKL model 
(X^ ed /dof= 1.1/1 04). Further XMM-Newton observations of this 
source are probably needed to resolve this issue. In Fig. [13] 
we show the spectrum of IGRJ08408-4503 accumulated over 
the entire exposure of the observation and fitted with the 
MKL+CUTOFFPL model. In this figure, we also show for com- 
parison the residuals obtained by fitting the same spectrum with 
a simple absorbed CUTOFFPL model (E cut =l 1 keV, fixed). The 
unfolded spectrum of IGR J08408-4503 is shown in Fig. HI 

We searched for pulsations in the power spectra of the XMM- 
Newton observation of IGRJ08408-4503 by using the same 
technique described in Sect. 14. II No significant (above 3cr level) 
signal was detected in these data. We determined an upper limit 
to the pulsed fraction of 25%, 30%, and 40% for periods in the 
range 0.3-50 s, 50-100 s, and 0.15-0.3 s, respectively (3a c.l.). 

5. Discussion 

We have presented the first deep pointed observations of the two 
prototypical SFXTs, XTE J 1739-302 and IGRJ08408-4503 in 
quiescence. 

The two sources exhibited a very complex timing and spec- 
tral variability, and we discuss them separately below. Here we 
also carry out a comparison between their quiescent and outburst 
emission. 

5.1. XTE J 1739-302 

The two XMM-Newton observations analyzed in Sect. 14. II show 
that the quiescent emission of XTE J 1739-302 is characterized 
by a number of low-intensity flares, taking place sporadically 
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IGRJ08408-4503 (Epic-PN) 




Energy (keV) 



Fig. 14. The unfolded spectrum of the XMM-Newton observation 
of IGRJ08408-4503. The best-fit model represented here is ob- 
tained by using an absorbed CUTOFFPL (cutoff energy fixed at 
1 1 keV) plus a MKL component. 



from a lower persistent emission level. The typical duration 
of these flares is a few thousands of seconds, and their X- 
ray flux is a factor of ^10-30 higher than the persistent qui- 
escent flux. During the time intervals when the source was at 
its lowest level of emission, we measured an X-ray flux of 
4.7 x 10 -13 erg/cm 2 /s (0.5-10 keV). This corresponds to a lumi- 
nosity of 4.1 x 10 32 erg/s at 2.7 kpc, and is among the lowest val- 
ues of X-ray luminosity measured from XTE J 173 '9 -302 (a factor 
of ^2 lower than the 3 a upper limit reported by Sakan o et alj 
2002 '), and compara ble with the lowest luminosity reported by 



I ... 

iRomano et alJ (l2009bl) ). The total dynamic range of the X-ray 
luminosity of XTE J 17 39-302 between outburst and quiescence 
isthus>10 4 (see Sect. IB- 

The hardness intensity diagrams and rate resolved analysis 
carried out in Sect. [4] showed that the variations in the X-ray 
flux measured during the XMM-Newton observations were ac- 
companied by a change in the spectral properties of the source, 
the source hardness ratio increasing significantly with the count 
rate. A fit to th e rate-resolved spe ctra with a CUTOFFPL model 
(-Ecut=F3 keV lSidoli et alll2009h implies that this behavior orig- 
inates from a change in the CUTOFFPL photon index, T, rather 
than in a variation of the absorption column density. 

These results indicate that the timing and spectral variability 
of XTE J 17 39-302 during the quiescence is qualitatively similar 
to that observed during the outbursts (see Sect. Q}. A change in 
the PL photon inde x with the X-ra y flux of XTE J 17 39-302 was 
first noticed bv lSmith et alj (Q998). These authors analyzed sev- 
eral different flares caught by RXTE, INTEGRAL and ASCA, 
and showed that the photon index of the hard (2-10 keV) X- 
ray emission from XTEJ1739-302 changed from 0.8, when the 
source X-ray flux was 2.4x 10 -9 erg/s (2-10 keV), to 2.0, when 
the X-ray flux was 1.6x 10 -9 erg/s (2-10 keV). However, these 
flares were also characterized by a significant change in the ab- 
sorption column density (from 3 to 37 x 10 22 cm -2 ). Similar val- 
ues of the PL photon index and absorption column density were 
reported for the outbursts observed with Swift from this source 
(see Sect.Q]). In the long-term monitorin g of XTE J 1739-302 car- 
ried out with SwiY uRomano et all (2009b) identified three differ- 
ent states in the source X-ray flux (Sect.[T|), and showed that each 
of these states could be characterized by a different PL photon 



index (in the range 0.5-1.6) and a different absorption column 
density (from 0.3 to 3.6 xlO 22 cm -2 ). Considering all of these 
results, it seems unlikely that there exists a single monotonic 
variation in both the absorption column density and PL photon 
index with the source intensity across its quiescent and outburst 
state. 

When the X-ray spectrum of XTE J 1739-302 from the en- 
tire XMM-Newton observation is considered, the interpretation 
of the emission from this source becomes even more compli- 
cated. In particular, this spectrum cannot be fit successfully by 
using a simple absorbed PL or a BB model. We indeed showed 
in Sect. 14.11 that these models provided only a poor fit to the 
data, and a more refined model was required. Because of the 
relatively low X-ray flux of the source, different models could 
equivalently describe the data. From a statistical point of view, 
the CUTOFFPL with a variable cutoff energy would be prefer- 
able, as it provided a very good fit to the data of both OBS 1 and 
OBS2 and requires a lower number of free parameters with re- 
spect to the others models reported in Table [2] However, this 
model would require a cutoff energy much lower (~4 keV) 
than the value determined when the source was in outburst 
(~13 keV). This low value for the cutoff energy might not be un- 
likely for an X-r ay pulsar (see e.g., t h e cases of X-Persei an d RX 
JO 146.9+61 2 fllHaberl et alJ (ll998h : lDiSalvo et al.l (fl998l and 
ICoburn et al.l (2002) for a review). However, we suggested that, 
given our relatively poor knowledge of the quiescent emission 
of XTE Jl 739-302 and of SFXTs in general, the applicability of 
other spectral models might be worth exploring. 

Besides the CUTOFFPL, the quiescent spectrum of 
XTE J 17 39-302 might require an additional spectral component 
at the lower energies (<2 keV). In Sect. 14.11 we showed that a 
fit to the spectra of OBS1 and OBS2 with a CUTOFFPL model 
and a fixed _E cut =13 keV would leave some evident structures 
in the residuals from the fit and the addition of a BB, or a MKL 
component can significantly improve the results. Even if these 
models require one more free parameter than the CUTOFFPL 
with a variable E cut , the results we obtained would then be in 
agreement with those found in the case of IGR J08408-4503. The 
spectrum of this source could, indeed, not be reproduced using a 
simple CUTOFFPL model (see below). 

A soft spectral component below 2 keV might be expected 
in the spectra of the SFXT sources, as this component is very 
common in bin aries hosting a NS a ccreting mass from a mas- 
sive companion. Iffickox et al.1(l2004l) showed that the detectabil- 
ity of this component depends mainly on the amount of absorp- 
tion in the direction of the sources. According to their results, 
in the most luminous objects (Lx— 10 38 erg/s) the soft compo- 
nent is produced by reprocessing hard X-rays from the NS by 
some optically thick material (e.g., an accretion disk), whereas 
for sources with lower luminosities (< 10 36 erg/s) the most likely 
origin of the soft component is the emission from either a pho- 
toionized or collisionally heated diffuse gas in the binary sys- 
tem or from the NS surface. In the case of XTE J 17 39-302, we 
found that BB emission from a relatively cold (~0.1 keV) and 
large (^100 km equivalent radius) region or from a much hot- 
ter (>1 keV) and less extended (<100 m) spot, or, alternatively, 
emission from an optically thin gas (MKL) provided equally 
good fits to the data. While the discussion above and the rapid 
variability observed in the SFXT would argue agains t the pres- 
ence of an accretion disk in these sources (see also. lBozzo et al.L 
2008), we also consider that emission from a small and hot spot 



Note that these authors used an absorbed power-law with exponen- 
tial high-energy cutoff (HIGHECUT*PL in XSPEC) to fit their data. 
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on the NS surface is inconsistent with models of accretion onto 
magnetic NS. The accretion flow in SFXT is, indeed, thought 
to be quasi-spherical, and expected to penetrate the NS mag- 
netosphere mainly by means of the Rayleigh-Tay l or an d the 
Kelvin-Helmholtz instability (see e.g., iBozzo et"all 120081 and 
references therein). In these circumstances, the size of the hot- 
spot over which the accretion takes place is expected to be in- 
versely proportional to the X-ray luminosity and might cover a 
substantial fractio n of the NS surface for Lx<l0 35 erg/s (see 
IWhite et al.L fl983. and references therein). 

Taking these results into account, emission from an op- 
tically thin and diffuse gas around the NS seems to be 
a more reasonable explanation of the soft spectral compo- 
nent of XTE Jl 739-302. According to this interpretation, the 
emitting region can be estimated from the normalization of 
the MKL component (see Table pj ) using the relation (see 
e.g., [Zurita-Heras & Walter! , 120091 and references therein): 

i?om=^/37VMKL/10- 14 (D/n H ) 2 ~10 13 aif cm, where AW 
is the normalization of the MKL component, D is the source dis- 
tance, nn~An/a, a is the binary separation, and ai3=a/10 13 cm. 
Therefore, the radius i? cm of the emitting region is compatible 
with the binary separation for a wide range of values of orbital 
periods similar to those measured in other SFXTs. We note also 
that the properties of this MKL component would be rather simi- 
lar to those derived from the spectru m of the SFXT AX J1845.0- 
0433 (IZurita-Heras & Waited 120091 

Another possibility that we investigated in Sect. 14.11 is the 
applicability of a model comprising a power law component and 
a partial covering to the X-ray spectrum of XTE Jl 739-302. We 
concluded that this model can also provide a reasonably good fit 
to the data. A similar model was proposed to interpret the qui- 
escent XMM-Newton spectrum of the SFXT IGR J 16207-5 129 
dTomsick et all 120091) and might provide support for clumpy 
wind in these sources. In this case, one would expect part of 
the radiation from the NS to escape absorption by the clumps 
local to the s ource and be affected only b y interstellar absorp- 
tion (see, e.g. [Walter & Zurita HerasLl2007l) . Finally, we showed 
that the COMPTT model also provided a reasonable fit to the 
XMM-Newton spectra of XTE J 17 39-302, and should thus be 
considered a valid alternative to the other models discussed 
above. We note, however, that neither the partial covering nor 
the COMPTT model could give an acceptable fit to the spectrum 
of IGR j '08408-4503 (see also below). Given the similarities be- 
tween the two sources, a spectral model that provides accept- 
able results in both cases should probably be favored (e.g, the 
CUTOFFPL+MKL model). 

5.2. IGRJ08408-4503 

The XMM-Newton observation of IGR J08408-4503, detected a 
very similar behavior to that discussed above for XTE J 17 39- 
302. In the light curve of IGRJ08408-4503, a number of rela- 
tively small flares were observed to take place sporadically on 
a timescale of few thousands of seconds and were characterized 
by an X-ray flux ^10-15 times higher than the underlying fainter 
persistent emission. The lowest X-ray flux that we measured 
from this source was ^3xl0~ 13 erg/cm 2 /s (0.5-10 keV) and 
corresponds to a luminosity of 3.3 x 10 32 e rg/s (at a distance of 
3 kpc) , comparable to the value reported bv lKennea & Campanal 
( 2006| . The total dynamic range in the X-ray luminosity of 
IGR J08408-4503 between outburst and quiescence is thus >10 4 
(see also Sect.Q]). 



As for XTE J 17 '39-302, the hardness intensity diagrams and 
the rate resolved analysis carried out in Sect. [4] showed that 
the variations in the X-ray flux measured during the XMM- 
Newton observation were also accompanied by a significant 
change in the spectral properties of the source. In contrast 
to the case of XTE J 1739-302, the rate-resolved spectra of 
IGR J08408-4503 could not be fit by using a simple CUTOFFPL 
model. We showed that an acceptable fit to the data could, in- 
stead, be obtained by introducing an additional relatively cold 
(fcT~0.08 keV) and extended (>100 km) BB component, or a 
MKL model (see Tables [3] and @). We note that the parameters 
measured for the MKL component are rather similar to those 
derived in the case of XTE J 17 39-302. An equivalently good fit 
was provided by the BMC model. This model has the same num- 
ber of free parameters as the CUTOFFPL+BB (see Tabled, and 
would predict similar properties for the temperature and the size 
in which the soft photons are produced (the normalization of the 
BMC model is defined as the ratio of the source l uminosity to 
the sq uare of the distance in units of 10 kpc, see e.g. lSidoli et all 
2009b, and references therein). Similar values of the fit param- 
eters were also obtained from the analysis of the spectrum of 
IGR J08408-4503 extracted by using the total available expo- 
sure time of the XMM-Newton observation (see Sect. 14.21 and 
TableSJ. 

According to the discussion in Sect. 15.11 a BB emission 
with these properties seems unlikely in the case of IGR J08408- 
4503 and thus we suggest that the CUTOFFPL+MKL model 
can provide a more reasonable description of the data. We note 
that in the XMM-Newton spectrum of the supergiant HMXB 
IGR J 16320-4751 a similar soft component was found that 
could be fit with a BB of 0.07 keV but w as attributed to a 
cloud surrounding the NS (Rodrigue z et all 1 20061) . Following 
the CUTOFFPL+MKL interpretation, the rate resolved analysis 
carried out for the observation of IGR J08408-4503 would indi- 
cate that the properties of the MKL component do not change 
significantly with the source count rate and the increase in the 
hardness ratio observed in Figs.QT]and[T0]is most likely due to a 
change in the CUTOFFPL photon index. Furthermore, no signif- 
icant variations in the absorption column density were revealed 
in the different rate-resolved spectra. This is similar to the results 
discussed above for XTE J 77 '39-302. 

As for XTE J 1739-302, a comparison between the results of 
the present XMM-Newton observation and the observations car- 
ried out in the same energy band (0.5-10 keV) when this source 
was in outburst does not indicate a clear correlation between 
the power law photon index, the absorption column density, 
and the source X-ray flux (see Sect.Q]). It is interesting that the 
soft component in this source detected by the XMM-Newton 
ob servation appears to have a different origin from that detected 
bv lSidoli et alT(l2009bh when IGRJ08408-4503 was in outburst 
(see Sect. Q]). On that occasion, the soft component appeared 
to be caused by thermal emission from a hot-spot on the NS 
surface. We note that, even if the soft component observed by 
XMM-Newton is interpreted in terms of a BB emission, the 
emitting region derived from the fit is considerably larger than 
the NS radius and it is thus unlikely that it is produced on the 
NS surface. 

Finally, for both XTE J 17 39-302 and IGRJ08408-4503, we 
investigated whether the harder spectral component detected in 
these sources might be produced by the X-ray emission from 
the NS supergiant companion. The time-averaged X-ray lu- 
minosity that we measured from these sources in quiescence 
matches quite well the luminosity expected from an isolated 
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OB supergiant or from collidi ng winds in a b i nary c ontaining 
OB supergiant stars (see e.g., iGudel & Naz el. 20091 for a re- 
view). However, this interpretation appears to be contrived for 
the following reasons. The X-ray spectrum of isolated or collid- 
ing wind binaries with OB supergiant stars is usually described 
well by a model comprising one or more thermal components 
(MKL in XSPEC, see e.g.. IGudel & NazelL l2009h . The softer 
MKL component has a typical temperature of ^0.2-0.7 keV, 
and is thus similar to those we detected in IGRJ08408-4503 
and XTEJ1739-302. This component is thought to be gener- 
ated by the shocks within the stellar wind. The hotter MKL 
component, extending up to several keV, can have a temper- 
ature as high as ~l-3 keV and is characteriz ed by a num- 
ber o f very prominent emission lines (see also iRaassen et al.L 
2008). This hard component is usually interpreted in terms of 
magnetically confined wind shocks, highly compressed wind 
shocks, or inverse Compton scattering of photospheric UV 
photons by relativistic p a rticles accelerated within the shocks 
(lAlbacete~ Colombo et all 120071) . Possible detections of a non- 
thermal X-ray emission from OB supergiant stars were re- 
ported only in only tw o cases, but they still lack confirmation 
(IGudel & Nazel I2009h . 

The X-ray spectra of XTE J 1739-302 and IGRJ08408-4503 
were reproduced well using a CUTOFFPL model, and no promi- 
nent emission lines were detected. The values of the photon in- 
dex, r, derived from XMM-Newton are also comparable to those 
obtained previously when the sources were in outburst, thus sug- 
gesting that a common mechanism produces their harder X-ray 
component. Furthermore, the relatively rapid X-ray variability 
(of period few thousands seconds) observed in the lightcurves 
of IGRJ08408-4503 and XTE J 1739-302 is not reminiscent of 
the typical X-ray variability of the OB stars, which, when 
present, takes place on longer timescales (tens of ks, see e.g, 
lAlbacete Colombo et all 120071) . 

We conclude that the harder X-ray emission from 
XTE J 1739-302 and IGRJ08408-4503 is most likely due to 
residual accretion taking place onto the NS at a much lower rate 
than during outburst. 



6. Conclusions 

The three XMM-Newton observations that we have analyzed in 
the present work, indicate that the quiescent spectra of the two 
prototypical SFXT XTE J 17 39-302 and IGR J08408-4503 are 
characterized by two different spectral components, one domi- 
nating the spectrum at the softer energies (<2 keV) and the other 
one being more prominent above 2 keV. 

The properties of the soft component (<2 ke V) could be rea- 
sonably well constrained in the case IGRJ08408-4503, where 
the absorption column density was relatively low (< 10 22 cm -2 ), 
whereas in the case of XTEJ1 739-302 the detection of this com- 
ponent is less significant. However, the similarity in the timing 
and spectral behavior observed in the quiescent state of the two 
sources argues in favor of adopting the same spectral model 
for both of them. We suggested that the model comprising a 
CUTOFFPL component at the higher energies plus a MKL com- 
ponent would provide a reasonable description of the data and a 
plausible physical explanation of the properties observed in the 
two sources (Sect.|5J- According to this interpretation, the MKL 
component would represent the contribution to the total X-ray 
emission of the shocks in the wind of the supergiant companion. 
The results of the fits with this model to the data of the three 
observations inferred a temperature of the MKL component and 



an emitting region comparable with the values found also in the 
case of the SFXT AX J1845.0-0433. 

Similar soft spectral components have been detected in many 
other HMXBs and SGXBs. In a few cases, the detection of 
a number of prominent emission lines in the high resolution 
X-ray spectra of these sources carried out with the gratings 
onbo ard Chandra and th e RGS onboard XMM-Newton (see 
e.g.,|\Vatanabe et all 120061) have convincingly demonstrated that 
these components are produced by the stellar wind around the 
NS, and proved to be a powerful diagnostic to probe the structure 
and composition of the stellar wind in these systems. The statis- 
tics of the present XMM-Newton observations is far too low to 
permit a similar in-depth study of the stellar wind in the case of 
XTEJ1739-302 and IGRJ08408-4503. Furthermore, because of 
the relatively low luminosity and the high absorption that charac- 
terize the emission of these sources in quiescence, observations 
at higher spectroscopic resolution are probably too challenging 
for the present generation of X-ray satellites, and the improved 
sensitivity of the X-ray spectrometers planned for future X-ray 
missions (e.g. IXO) is probably required to firmly establish the 
presence of a soft spectral component in the quiescent emission 
of the SFXT sources and shed light on its nature. 

If the harder X-ray emission (2-10 keV) detected from the 
XMM-Newton observations of XTE J 17 39-302 and IGRJ08408- 
4503 was really produced by residual accretion as we argued in 
the previous section, then the accretion process in these sources 
would take place over more than 4 orders of magnitude of X- 
ray luminositjQ This is similar to the resu lts reported for the 
SFXT IGR J 17544-26 19 dRampv et all [2009b and, possibly, for 
the SFXT SAX J1818.6-1703 (in the latter case the origi n of th e 
lowest quiescent emission remain s unclear, Bozzo et all 12009). 
In the case of IGR J 17544-26 19, |RamEy_etaL| (|2009|) ascribed 
the high dynamic range in the X-ray luminosity to the accretion 
of clumps from the wind of the supergiant star with a high den- 
sity contrast with respect to the surrounding homogeneous wind. 
However, it was also suggested that a similar variability might 
result fr om the transition ac ross different accretion regimes onto 
the NS dBozzo et all 12008b . We note that a similar scenario can 
be envisaged for interpreting the variations in the X-ray flux ob- 
served during the multiple small flares detected in the present 
observations. Even though they took place at a much lower lumi- 
nosity level than the brightest outbursts (a factor of ~10 3 -10 4 ), 
our analysis showed that all these events shared a number of sim- 
ilar timing and spectral properties. In particular, the timescales 
on which the smaller flares develop is comparable with the de- 
cay timescale of the source luminosity during the outbursts (see 
Sect. [TJ, and the spectral photon indices and absorption column 
densities measured from the XMM-Newton observations are also 
in qualitative agreement with those reported previously when the 
sources were observed at a much higher X-ray luminosity level 
(see Sect. Q]). It is thus most likely that the transitions between 
the lower quiescent states and the small flares detected by XMM- 
Newton from XTE J 17 39-302 and IGRJ08408-4503 might have 
been triggered by the same mechanism that sometimes gives 
rise to the brightest outbursts (i.e., the accretion of clumps from 
the stellar wind and/or the transition between different accretion 
regimes of the NS, see Sect.[T|i. 

In contrast to the case for the SFXT IGR J18483-031 1, we 
did not detect any pulsation in the quiescent emission of either 
XTE J 1739-302 or IGRJ08408-4503, and provided in Sects. |4~T1 



That we did not find any evidence for X-ray eclipses in the three 
XMM-Newton observations that we analyzed (this is unlike the case of 
IGR J 16479-45 14, see Sect.fTJ is also consistent with this interpretation. 
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and 14.21 the corresponding upper limits to the spin periods and 
pulsed fractions we were able to infer from the present data. 

Deep pointed observations of SFXTs in quiescence are still 
required in order to understand the origin of the peculiar X- 
ray variability of these sources and distinguish between different 
models proposed to interpret their behavior. 
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